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We study the nuclear modification of the scalar QCD susceptibility, calculated as 
the derivative of the quark condensate with respect to the quark mass. We show that 
it has two origins. One is the low lying nuclear excitations. At normal nuclear density 
this part is constrained by the nuclear incompressibility. The other part arises from 
the individual nucleon response and it is dominated by the pion cloud contribution. 
Numerically the first contribution dominates. The resulting increase in magnitude 
of the scalar susceptibility at normal density is such that it becomes close to the 
pseudoscalar susceptibility while it is quite different in the vacuum. We interpret it 
as a consequence of chiral symmetry restoration in nuclei. 
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I. INTRODUCTION 

The chiral phase transition is a hot topic of QCD. The attention has focused in particular 
on a baryonic rich environement, the simplest one being ordinary nuclei. In view of the 
difficulties of the lattice methods to treat baryonic matter, studies based on models have 
revealed their interest. In the first stage the order parameter, which is the quark condensate, 
was investigated [l|, Q|. With the realization that the amount of restoration is large, as the 
order parameter has decreased by about 35% at ordinary nuclear density, the interest has also 
focused on precursor effects linked to the partial restoration of chira 
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symmetry in the form 
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attract attention is the question of the susceptibilities in QCD[8, Q|. In the broken symmetry 
phase, where the order parameter introduces a privileged direction, the susceptibility is 
splitted as in magnetism in the parallel susceptibility, along the magnetization axis, and the 
perpendicular one. In QCD these are the scalar (xs) and pseudoscalar (xps) susceptibilities 



related to the fluctuations of the scalar and pseudoscalar quark density, respectively. 

The in-medium QCD susceptibilities have been discussed by Chanfray and Ericsonj^]. For 
the scalar one they used the linear sigma model which provides the coupling of the quark 
scalar density fluctuations to the nucleonic ones through sigma exchange. They ignored 
the coupling to the pion density fluctuations, expected to have a smaller influence. Their 
approach is basically a dispersive one, with the introduction of the in-medium scalar spectral 
function. In terms of graphs their effect corresponds to the one of Figured and the dressing 
by the pion lines to that of Figure [2j 

The present work uses a totally different approach which relies on the very definition of 
the longitudinal susceptibility as the derivative of the order parameter with respect to the 
perturbation responsible for the explicit symmetry breaking. In magnetism it is the external 
magnetic field. In QCD it is the quark mass and we have the generic definition of the scalar 
susceptibility: 

Xs = ^ q <(iq> > (1) 
where q is the quark field. The two methods are known to be equivalent in their principle. 
However in practice, in the dispersive approach, truncations are made in the intermediate 
accessible states. In fact our new derivation incorporates terms as well as interference effects 
which were previously absent. The present study is also free from the specific features of 
the linear sigma model. 

In Section UH we show how Eq'JT]) leads to a natural decomposition of xs into a vacuum 
contribution xs{ vac -), a contribution noted Xs whis is related to the nucleonic excitations, 
and a contribution noted ^ uclear which is related to the nuclear excitations. In Section 
IIIII we use the Fermi gaz model to estimate ^ uclear and we compare the result with the 
one obtained previously with the linear sigma modelQ. We conclude that ^ uclear depends 
essentially on the zero momentum particle-hole propagator which we relate to the nuclear 
incompressibility. In Section 4 we study the nucleonic contribution Xs an d we find that it 
is dominated by the pion cloud, the quark core giving only about 6% of Xs- The model 
dependance of the pion cloud contribution is moderated by the fact that, within a factor of 
order two, its value is fixed by the leading non analytic piece of the sigma term. In Section 
El we present some numerical estimates based on the results of Sections 2, 3 and we discuss 
their implications concerning the restoration of Chiral symmetry in nuclei. 



II. NUCLEAR SUSCEPTIBILITY 



Since our aim is to evaluate the modification of the susceptibility with respect to its 
vacuum value we note < qq(p) > the in-medium value of the quark condensate and we 
introduce: 

XS(P) = ^ < qq( P ) > (2) 

where p is the nuclear density. 

For a collection of independent nucleons, the in-medium condensate < qq(p) > writes: 



< qq(p) >=< qq(vac.) > 



j ^nPs 
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Here ps is the nucleon scalar density, < qq(vac.) > denotes the vacuum expectation value 
of the condensate and T, N is the nucleon sigma term: 

S^v =< N\ [Q s ,Qs] \N >= 2m q J dx < N\qq(x) - qq(vac.)\N > . (4) 
Using the Gellmann, Oakes and Renner relation: 

fl m l = ~ 2m q < qqivac.) >, (5) 
the in-medium condensate expression (jHJ) can also be written in the form: 

< qq(p) >=< qqivac.) > +— . (6) 

The result above follows from the Feynman-Hellman theorem which relates the condensate 
to the thermodynamical grand potential per unit volume lo = e — pp through: 

where the derivative has to be taken at constant baryonic chemical potential p (which 
controls the density p). As an illustration in a free Fermi gas, one has after substracting the 
vacuum energy: 
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with E p = y/p 2 + M 2 . The medium contribution to the condensate is obtained as: 
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< qq(p) >-< qq{vac.) >= - — - — = p s — , (9) 
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where the scalar density is defined by: 



^ = 4 /(^f e( "-^ (10) 



while the ordinary density is: 



P 



with p 2 F = p 2 - M 2 . 

Note that in Eq.(PJ the contribution of the delta function from the derivative of the 
Heaviside function vanishes. It will not be the case for the susceptibility. Starting from 
Eq.fjnj) one gets: 

/ \ 9 _ d {E N \ H N dps 

Xsip) = d^ q < qq{vac) > +Ps d^ q ) + 2^ 0< (12) 

which contains three contributions: 

1. The derivative of < qq(vac. >, which is the vacuum susceptibility xs(vac). Its evalu- 
ation would require a non perturbative QCD model, which is outside the scope of this 
work. So we focus on the difference Xs{p) ~ Xs( vac -)- 



2. The derivative of S^r /2rn g , which is in fact the nucleon scalar susceptibility Xs- 
This follows from the relation between the sigma term and the condensate: 

= 2m q J dx < N\qq(x) — qq{vac.)\N > . (13) 



Thus: 

Xs 



N - 9 ■ f dx<N\qq{x)-qq{vac.)\N >=^-[^-\. (14) 
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Therefore this second term, which writes psXs i can be interpreted as the individual 
nucleon contribution to the nuclear susceptibility. 



3. The derivative of ps which gives the third term, noted ^ uclear : 

^nuclear _ 9p$ . , 

We shall see that it represents the effect of the nuclear excitations by contrast to the 
second term which is due to the nucleon excitations. 

In summary we are going to study separately the two terms of the quantity 

Xs( P ) ~ Xs(vac.) = PsX % + Xs udear - (16) 

III. NUCLEAR EXCITATION CONTRIBUTION 

We first examine the term y^ udear in Eq. (|16jl . We use the free Fermi gas, for which 
the scalar nucleon density is given in Eq. lfTUl) . Taking its derivative with respect to the 
quark mass at constant chemical potential we get two terms, the second one arising from 
the derivative of the Heaviside function: 
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The first term in Eq. (|T7jl represents the polarization through nucleon-antinucleon (NN) 
excitations of a relativistic Fermi gaz submitted to a scalar perturbation (see Appendix). 
As it vanishes in the non relativistic limit we shall neglect this NN contribution. 
From the second term in Eq. (fT7|) we obtain: 

In the non relativistic limit (/i ~ M) the parenthesis reduces to (—2Mp F /7i 2 ) which is 
actually the particle-hole polarization propagator H p h(q) of the non relativistic free Fermi 
gaz taken in the static situation, i.e., for q = 0, and taking the limit of vanishing three 
momentum (q — > 0). We have then: 



The presence of U p h indicates that the origin of this term lies in the nuclear excitations. 



A. Comparison with y™ Lclear obtained in the sigma model 

In the sigma model the scalar susceptibility is related to the a propagator jEj. In the 
medium this propagator is modified by the particle-hole insertions, which gave: 

X7 dear = 2 <9 ' 9( ? 2 C ° >2 4 n *(°> 0) (2°) 
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where II 5 is the full scalar particle-hole polarisation propagator and g a is the sigma nucleon 
coupling constant. In order to link the two expressions of ^ uclear we first need to evaluate 
the nucleon sigma term £jv within the linear sigma model. It is built of two pieces. The 
first one corresponds to the exchange of a sigma between the condensate and the nucleon as 
illustrated on Figure Ek- The second one is the contribution of the pion cloud of the nucleon 
which comes into play as two pions exchange, as shown in Figure |2]b. The result is 



2 

,.v = + = /,m^ + ^ dx < iV|0 2 (f)|iV > . (21) 



In the context of our previous workQ], the pionic contribution Y? N did not appear naturally 
and thus was ignored. To make the comparison meaningful we should then retain only the 
sigma exchange part and insert it in our expression (fT9|) which leads to: 

x7 dear (°) = ^4iWo, o) = 2 <^-)> 2 ^ np/t(0) 5) (22) 

This is essentially the result of the sigma model calculation (see Eql2"0|). provided we replace 
the full n s by the free Fermi gas expression Yi ph . The present derivation is more satisfactory 
in the sense that it does not rely on the sigma model to derive the coupling between the 
quark density fluctuations and the nucleon ones. In particular it incorporates other couplings 
than through sigma exchange, such as the two pion exchange term of T>n- Moreover the 
interferences between the various components of S^r are automatically incorporated in the 
crossed terms of T? N . One of these interferences, that is the S^E^ term, is illustrated on 
Figure |U 

On the other hand talking also about the limitations of the present work, it applies to a 



free Fermi gas (its generalization is in progress) , while our previous approach did not restrict 
to this situation. In the dressing of the sigma line by ph states, the full ph propagator 
entered. The latter is to some extent constrained by nuclear phenomenology [a| and the 
argument goes as follows. At normal nuclear density there is no distinction between the 
scalar and ordinary density operators. Now, at ordinary density, the ph progator U p is the 
response of the system to a perturbation which couples to the nucleon density. In other 
words it is the nuclear compressibility, with the relation: 

n P = (23) 

where the nuclear incompressibility K is related to the energy per particle E/A by 

d ( 9 d E\ 

Its experimental value at the saturation density p is m the range 200 -r300 Me V This 
value is compatible with the free Fermi gas value computed at the same density. Even though 
this agreement may result from a cancellation between several many-body influences, such 
as the effective mass and residual interaction effects, it justifies to some extent the use of 
the free Fermi gas model. 



IV. NUCLEON SCALAR SUSCEPTIBILITY Xs 

For a structureless nucleon one has of course Xs = but we do not make such a restriction. 
In general Xs because the true nucleon responds to a scalar perturbation by adjusting 
its internal structure. One can estimate this response using models, for instance the MIT 
bag model, but it could also be extracted from lattice calculations when the latter are done 
at realistic quark mass values. 



A. Valence quark contribution to Xs 

The scalar susceptibility of a free nucleon has been introduced by Guichonjllj] in another 
context, concerning a pure nuclear physics problem, that is the question of saturation of 
nuclear matter. In his quark-meson coupling model the saturation follows from the response 



of the nucleoli to a scalar field. In the bag model the total scalar charge of the bag defined 
as: 

Qs = J dx<N\qq(x)\N>, (25) 
depends on the quark mass. The linear term in the quark mass expansion of Qs' 

Qs(m g ) = Q s (0) + xs {Bag)m q + ■■■ (26) 

defines the susceptibility of the bag. It is found to be Xs {Bag) ~ 0.5-R ~ 2.5 10 -3 MeV -1 , 
where the numerical value corresponds to a bag radius R = 0.8/m. It turns out that this 
susceptibility due to the quark structure can stabilize the scalar nuclear field and provide 
a mechanism for saturation[3|. In our case however its contribution to the nuclear suscep- 
tibility is negligible (see Section EJ) with respect to the one due to the nuclear excitations. 
This is not a surprise since the nucleonic excitations, which control Xs(Bag), are much 
higher than the nuclear ones. What is interesting however is the positive sign of this term. 
Its origin is rather obvious: when its mass increases the quark becomes less relativistic. This 
tends to increase the scalar charge of the nucleon as the quark scalar density qq involves the 
difference of the large and small components of the quark spinor. 

Another important point is that Xs(Bag) only includes the valence quark contribution 
and not those from the pion or sigma clouds. The argument about the energy of the corre- 
sponding excitations may not apply to the pion cloud contribution in view of the small pion 
mass. It is therefore interesting to evaluate the corresponding susceptibility. 



B. P ionic contribution to Xs 

The nucleon sigma commutator is largely influenced by the presence of the pion cloud. 
If the nucleon remains unexcited after pion emission and in the heavy baryon limit the 
corresponding contribution is equal to (3| 

where <j) is the pion field and F(q) the tcNN form factor. With this explicit expression of 
S^r it is straightforward to calculate the derivatives with respect to m q (or m%) involved in 



the corresponding susceptibility: 
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For a monopole form factor F(q) = A 2 / (A 2 + q 2 ) we find: 

_9m| <TO( B ac.)> 2 ^\ 2 / A V = _ 4 10 -2 Me y-, 



tf W = ( — j = -4HTW (29) 

where the numerical value corresponds to A = 5m„-. This value is about 15 times larger 
than the susceptibility due to the quark bag structure estimated in the previous section. 

A rough estimate of the pionic contribution can be obtained if the form factor F(q) is 
omitted in Eq.JSHJ), which corresponds to the limit A — > oo. In this case Xs i 71 ) can be written 
in terms of the leading non analytical term of the sigma term S^" 4 according to: 

N 2 < qqiyac.) > 2 d Y, L N NA < qq{vac.) > 2 LNA 

XsW = Ti *4^T = TiK E " ' (30) 



where 
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= -— 1^1 mi * -23MeV. (31) 



This approximation leads to Xs( 7r ) — — 8.6 10~ 2 Me^ _1 , that is about the double of the 
value obtained with the form factor. 

There is also a contribution to the nucleon sigma term due to intermediate ttA states 
(see Figure El b). In the narrow width approximation it reads: 

tA 3ml fg A \ 2 4 f g wNA V [<*>,( q 4 q 4 



S5r -i6*H/J qUvJ Jo *[h 2 k+v + kk+a 9 )J f(9) (32) 



with oj q = yj q 2 + m 2 and A q = Ma — M N + q 2 /2M&. The corresponding contribution to the 
susceptibility, Xs ( 7r ^)i * s obtained from the derivative of the above expression with respect 
to m 2 . The presence of the large energy denominator Ma — makes it less sensitive 
to the pion mass. Numerically with the ratio tva / Qtt n n = v/72/5, we find Xs — 



—1.410 2 MeV The overall susceptibility of a free nucleon due to the pion cloud is thus 
about -5.4 lO^MeV- 1 . 

Our conclusion on the nucleon scalar polarisability is that, as the electric one, it is domi- 
nated by the pion cloud. Within a factor of two it can be expressed in a model independent 
way in terms of the leading non analytic part of the nucleon sigma term. In this context 
it is legitimate to wonder why the pionic susceptibility Xs M which is so dominant in this 
problem does not also dominate the nuclear saturation problem where instead it is totally 
ignored. The answer lies in the chiral properties of the scalar field responsible for the nuclear 
attraction that we have studied in 14]. We have stressed that this field has to be chiral 
invariant rather than a, the chiral partner of the pion. It couples derivatively to the pion. 
Therefore, in the long wave length and static limit, the pion cloud is weakly coupled to this 
nuclear scalar field. 



C. Interpretation of the pion cloud contribution in the sigma model. 

It is interesting to look at the nuclear susceptibility PsXsi 71 ) i n the framework of the 
sigma model. The sigma, which transmits the quark fluctuations, is dressed not only by the 
particle-hole excitations but also by the two pions excitations, as shown in Figure El Since 
we are interested in the modification of the susceptibility with respect to the its vacuum 
value, at least one of the pions in this graph has to be a nuclear one. So it is dressed by 
particle-hole insertions. To lowest order in the density the contribution of this graph to the 
nuclear susceptibility is: 

f%™i ±f%\f*J J (2tt) 4 (q 2 - mlf 

Here the subscript L in Ul indicates the spin longitudinal character. In the static limit the 
domain where II(g , q) has a non vanishing value is pushed to zero energy. In this case the 
pions do not carry energy so the integral over q° which then involves only Il(g , q) reduces 
to: 

d *>TT / fdq f ( 1u)\ ImU L (ujq) 



7T J q 2 — L) 2 + VlTj 



J du (--ImIl L (uj q)^j = p J duR L (uq) = pS L (q), (34) 



where Rl{uj q) = —ImXl L (uj q) /tt p is the nuclear longitudinal spin-isospin response and Sl 
its integral over energy. For a free Fermi gas one has: 



S L (q)=e(q-2p F ) + e(2p F -q) 
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Ignoring the Pauli blocking effect the quantity Si(q) reduces to unity, which leads to: 
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(36) 



This is exactly the nucleonic polarisability arising from the pion cloud multiplied by the 
density. In fact the mere comparison of the many-body graph of Figure El with the one 
of Figure El which represents the free nucleon susceptibility leads to this conclusion. To 
leading order the evaluation of the nuclear QCD scalar susceptibility arising from the 2tt 
intermediate states does not need any calculation as it is simply related to the nucleon one, 
if Pauli blocking is ignored. 



V. NUMERICAL ESTIMATES 



We have now all ingredients to proceed to the numerical evaluation of the in-medium 
modification of the scalar susceptibility. To fix the idea we shall use IV — 45 MeV and 
K ~ 230MeV~ 1 . At normal density the contribution of the nuclear excitations to the 
susceptibility is then: 

Xs udear (po) = -8.2 10 5 MeV 2 . 

Turning now to the nucleonic participation to the scalar susceptibility we have to take into 
account the Pauli blocking which reduces the pionic cloud contribution from irN intermediate 
states by about 25% at normal density, bringing the in-medium nucleon susceptibility to the 
effective value Xs — — 4.9 10~ 2 MeV _1 . We multiply it by the nuclear density which gives 
PoXs — —6.8 10 4 MeV 2 . This number is smaller than the nuclear excitation contribution. 
Altogether the scalar susceptibility of nuclear matter at normal density is: 



Xs(po) = Xs{vac.) - 8.910 5 MeV 2 . 



(37) 



It is interesting to give a scale to compare this nuclear modification of the scalar susceptibil- 
ity. The susceptibility of the vacuum xs( vac -) is n °t known but due to the large mass of the 
scalar meson, it is certainely much smaller than the pseudoscalar susceptibility xps(vcic). 
The latter is actually dominated by the Goldstone boson, i.e. the pion, which allows its 
evaluation. Chanfray and Ericson give the following expression Q]: 

XPsivac) = - ^"f^ ^ - "1-3 10 6 MeV- (38) 
From this we infer that: 

• Xs(vac.) can reasonnably be neglected on the right hand side of Eq. ljHTj) . 

• at p = po the nuclear scalar susceptibility is comparable to the vacuum pseudoscalar 
susceptibility XPsiyac). 

Moreover Chanfray and Ericson have shown that Xps(p) follows the density evolution of 
the quark condensate, i.e., at normal density it has decreased by 35%, which brings it to 
Xps(Po) — -910 5 MeV 2 . This is nearly equal to the value Xs(Po) — -8.9 10 5 MeV 2 which 
we get when we neglect Xs(vac) in Eq. pTj) . This implies that the scalar and pseudoscalar 
susceptibilities, which are so different in the vacuum, become nearly equal at normal nuclear 
density, a feature normally expected only near the chiral phase transition. As our evaluation 
uses the value of the free nucleon sigma term, this convergence of Xs(po) and Xps(Po) toward 
a common value must be taken with a grain of salt due to the possible medium renormal- 
isation: S^v — > Sjv(p ), which we have not taken into account in this work. However this 
manifestation of the restoration of chiral symmetry is so spectacular that we do not expect 
it to be totally destroyed by these renormalisation effects. A systematic investigation of this 
problem, as well as the extension of this study at larger densities, deserves further work. 

In summary we have found that the two QCD susceptibilities, namely the scalar and 
pseudoscalar ones, undergo a strong modification in the nuclear medium in such a way they 
become close to each other already at normal nuclear matter density. It is a spectacular 
consequence of partial chiral restoration which may have consequences in processes involving 
higher densities. 



Appendix 



Here we want to interpret the first term in Eq. (|17|) in terms of nucleon antinucleon 
excitations. For this we can consider a relativistic free fermi gas described by a hamiltonian 
H and add a pertiirbation|lfil| 

\W = xj drtpip = xj dfps(f) (39) 

which changes the nucleon mass by the amount 5M = X. If we note |A > the ground state 
of the system in presence of the perturbation then, by the Feynamn Hellman theorem, we 
have 



^ < A| J df Ps (r)\X > ^ = V^< X\p s (0)\X >| A=0 = 2 < A = Q\W Eq ^ h W\X = >, 

(40) 

where we assume < A|A >= 1 for simplicity and V is the volume of the gas. To simplify we 
note 

^<A|p 5 (0)|A>| A=0 = ^. (41) 

From the cannonical field expansion, the part of W which produces the NN intermediate 
states is 

W N n = J dp^[u{p)v{-p)b ] {p)d ] (-p)} + h.c. 

= - f dp?-f[tf($)d\-p)] + h.c. (42) 
J h p 

where 6, d are respectively the destruction operators of the nucleon and antinucleon. So the 
NN contribution to the R.H.S' of Eq.fgOI) is: 
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where we have used 



< A = 0\tf(p)b(p)\ A = >= 46(0)9( PF - p) 
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(44) 



The infinite term is independent of the density so it drops out when we substract the 
vacuum contribution. Since the perturbation (|3H is equivalent to a change 5M = A of the 
nucleon mass, we can write 
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which is the first term in Eq. lfTTjl . In other terms the derivative of the scalar density at fixed 
density is entirely due to the NN excitations. 

Note that this does not allow us to say that the term 



Sat dps 



2m q dm q 
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in Eq.tJTnj is the nuclear susceptibility due to iVTV excitations. To reach such a conclusion 
we should start with a perturbation of the form 



m q J df[uu(r) + dd{r)\ = m q W, 



(47) 



that is the true mass term of QCD. Defining the nuclear susceptibility (per unit volume) as 



nuclear 
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2V dm~ <rrLq][ J + dd(r)}\m q > 



(48) 



m„=0 



where \m q > denotes the nuclear ground state in presence of the quark mass term p7l) then 
the Feynman-Hellman theorem gives 



nuclear 
AS 



2 < m q = QlW^-^-Wlm, = > . 



(49) 



To compute this quantity using the Fermi gas approximation we need to know the matrix 
element < iVlW^iV > for the ph excitations and < NN\W\vac. > for the NN excitations. 
There is no problem with the first one since we know it from the nucleon sigma term Y, N . 
On the other hand the second one is essentially unknown. We can quantify this problem by 
introducing the NN and NN scalar form factors in the standard way: 

m q < N(p')\uu(0) + dd{0)\N{p) > = S NN [{p-p') 2 ]u{p')u{p), 
m q < N(p')N(p)\uu(0) + dd(0)\vac. > = S N *[(p + p') 2 }u(p')v(p) (50) 

From the definition of the nucleon sigma term: 



< N(0)\N(0) > 



< N(0)\ / dfm q [uu{f) + dd{r)]\N(0) >, 



(51) 



we get S (0) = Sjv and by the crossing rule S and S are the same function. So we 



can write 



< N(p')N(p)\uu(0) + dd(0)\vac. >= —g[(p + p') 2 ]u(p')u(p) 



(52) 



where we have defined g(x) = S (x)/S (0). A straightforward calculation then leads to 



the following expression[l' 



for the NN contribution to the nuclear susceptibility: 



Xs udear (NN) 



Y 2 

N 



2m 2 (2tt) 



PF 



dp g{AEl 



2 p z 



g(AM 2 



2 Sjv dps 



2m q dm q 



(53) 



NN 



which differs from (jjQ| by the factor |^(4M 2 )| . This factor is probably very small because 
the transition vacuum — ► NN through the one body quark operator uu + dd is suppressed 
as compared to the elastic transition iV — > N. 
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Figure 1: Modification of the a propagator by the particle-hole polarisation propagator. The cross 
represents the condensate 




Figure 2: Modification of the a propagator by the in medium 2ir propagator. The cross represents 
the condensate 




Figure 4: Interference between a and 27r exchange. 
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Figure 5: Pionic contribution to the nucleon susceptibility. 




Figure 3: Nucleon sigma term in the a model 



